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^resonances. We compute the next-to-leading order QCD corrections to graviton plus jet hadro- 
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^^The QCD corrections are sizable and lead to a significant reduction of the scale dependence. We 
^_^p]resent numerical results for cross sections and distributions, and discuss the uncertainty from 
^nparton distribution functions and the ultraviolet sensitivity of the theoretical prediction. 
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1. Introduction 

The search for new physics at the TeV-scale is one of the major tasks for current and future 
high-energy physics experiments. Models with extra space dimensions and TeV-scale gravity 
address the problem of the large hierarchy between the electroweak and Planck scales, and 
predict exciting signatures of new physics that can be probed at colliders jy]. 

In the D = 4 + 5 dimensional model proposed by Arkani-Hamed, Dimopoulos and Dvali 
(ADD) 0, the standard model (SM) particles are constrained to a 3 + 1 dimensional brane, 
while gravity can propagate in a 4 + 5 dimensional space-time. For simplicity, the additional 
(5-dimensional space is assumed to be a torus with common compactification radius R. In such 
a model, the 4-dimensional effective Planck scale Mp is related to the fundamental scale Mg 
by §: 

Ml = 8ttR s M^ +2 . (1.1) 

For a large compactification radius R it is thus possible that the fundamental scale is near 
the weak scale, Ms ~ TeV. In the ADD model, deviations from the standard Newton law 
of gravity are predicted at distances around R ~ 0.83 x io~ 16 + 30 / 5 mm x (2.4TeV/M s ) 1+2 / 5 . 
Current terrestrial test of gravity exclude R > 37(44) /im for 5 = 2(1) |3L which, using Eq. (|1.1|), 
translates into Mg > 3.6 TeV for 5 = 2. Further constraints have been derived from astrophysics 
and cosmology, in particular for 5 < 4, but they can be evaded in specific models @ and do not 
lessen the importance of collider searches for extra dimensions. 

The D = 4 + 5 dimensional graviton corresponds to a tower of massive Kaluza-Klein (KK) 
modes in 4 dimensions. The interaction of these spin-2 KK gravitons with SM matter can be 
described by an effective theory H |], [jj with the Lagrangian 

£ ^ = -^E G ^ )T ^' (1-2) 

n 
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where the massive gravitons are labeled by a (^-dimensional vector of integers, n = (ni, --,ns), 
Mp = Mp/y/8x ~ 2.4 x 10 18 GeV is the reduced 4-dimensional Planck scale, and T^ v is the 
energy-momentum tensor of the SM fields. The Feynman rules that follow from Eq. ( |1.2|) can be 
found in Refs. || [/J. The individual KK resonances have masses equal to mr^ = \n\/R and thus 
the mass gap between neighboring modes Am = R^ 1 is small for 5 not too large. Quantitatively 
one finds Am « 20 keV, 7 MeV and 0.1 GeV for Mg = 1 TeV and S =4, 6 and 8, respectively §. 
The discrete mass spectrum can thus be approximated by a continuum with a density of states 
diV = p{m)dm ]|, ||], where 

, > M P * , , _ 2vr <5 / 2 , , 

p(m) = Ss-ij^m , and S 5 _i = r ^^ 2 y ( L3 ) 

Inclusive collider cross sections, where one sums over all accessible KK modes, are obtained from 
a convolution of the cross section for an individual KK mode of mass to, da m , with the mass 
density function p(m) Q1.3Q, da /dm = p(m)da m . Although each individual graviton couples 



2 

to SM matter with only gravitational strength oc 1/Mp, see Eq. (11.21), and thus da m oc 1/M P , 

— 2 

inclusive collider processes are enhanced by the enormous number of accessible KK states oc M P 
(|l.3|). The factors Mp cancel in da /dm = p(m)da m , leaving an overall suppression of only 
Mg . If the fundamental scale Mg is near the TeV-scale, graviton production can thus be 
probed at present and future high-energy colliders. 

Both virtual graviton exchange between SM particles and real graviton emission provide 
viable signatures of large extra dimensions at colliders. Since the coupling of gravitons with 
matter is suppressed oc 1/Mp, direct graviton production gives rise to missing energy signals. 
Searches for graviton production have been performed in the processes e + e~ — > 'y(Z) -\- E miss at 
LEP and pp — > 7(jet) + p™ lss at the Tevatron. The combined LEP limits || read Mg > 1.60, 
1.20, 0.94, 0.77, 0.66 TeV, for 8 = 2,- • • ,6 respectively, while Tevatron searches exclude Mg > 
1.40, 1.15, 1.04, 0.98, 0.94 TeV, for S = 2,- • • ,6 respectively |, 0. Searches for the process 
pp —> jet + p™ lss at the LHC will be able to extend the sensitivity to the fundamental scale Mg 
into the multi-TeV region [0, 00,^]. 

Current analyses of graviton production at hadron colliders are based on LO cross sections, 
which are subject to large theoretical uncertainties from the choice of renormalization and fac- 
torization scales. In this paper we present the first calculation of the NLO QCD corrections to 
graviton production in the process pp/pp — > jet + G at hadron colliders, results for the QCD 
corrections in the photon channel have been presented in [15]. The NLO cross sections lead to 
significantly more accurate theoretical signal predictions and thereby more accurate constraints 
on Mg or, in the case of discovery, will allow to probe the model parameters. The NLO calcula- 
tion also enables us to properly study the relative importance of multi-jet final states, which has 
been addressed in a recent calculation of the (tree-level) graviton plus di-jet cross section [16]. 
We note that NLO QCD corrections to the hadro-production of lepton and boson pairs in models 
with large extra dimensions have been presented in a series of recent papers [17, 18]. 

The remainder of this paper is structured as follows: In section |2| we present some details of 
the NLO QCD cross section calculation for pp/pp — > jet + G. Numerical results for the Tevatron 
and the LHC are presented in section ||. We conclude in section [|. More details of the calculation 
and selected formulae are presented in the appendix. 
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2. Calculation 



The LO cross section for graviton plus jet production receives contributions from the partonic 
processes 

qq -> gG, qg -> qG and gg ^ gG . (2.1) 

The relevant Feynman diagrams are depicted in Fig. |]. The LO partonic cross sections of the 
processes (|2.1| ) have first been presented in Ref. ||; the corresponding helicity amplitudes are 
discussed in the appendix. 

// ~ 
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Figure 1: The LO Feynman diagrams for gg,qq — > G+jet. The gq channel is related to the qq channel 
by crossing and not shown. 



The NLO cross section consists of virtual corrections, real-emission contributions and a 
collinear term, which is a finite remainder of the factorization of collinear singularities into the 
parton distribution functions (PDFs). We use dimensional regularization in d = 4 — 2e 
dimensions to regulate the ultraviolet (UV) and infrared (IR) divergences, and apply the dipole 
subtraction scheme [pi] ] to cancel the infrared singularities. The UV divergences are removed 
by renormalization of the QCD coupling a s in the MS-scheme. 

We have performed two independent calculations of the virtual corrections, described in 
more detail below, and have checked gauge invariance and Ward identities arising from general 
coordinate invariance, see Ref. [|l6| for more details. The numerical implementation of the real- 
emission contributions is based on MadGraph ]2^| and MadDipole Some details of the NLO 
calculation are provided below. 



2.1 Virtual corrections 

The virtual corrections to pp/pp — > jet + G arise from the interference of the Born and one- loop 
amplitudes. Example one-loop Feynman diagrams are depicted in Fig. ||. The Feynman rules 
for the graviton interaction with the SM fields can be found in Ref. ]|, |?], [T^] 1 . For external 
gluons we choose the light-cone gauge to avoid introducing external ghost lines, as in the HELAS 
convention [24]. The internal gluon propagators are evaluated in the 't Hooft-Feynman gauge, 



1 Note that there are different sign conventions in the definition of the covariant derivative, which leads to a 
sign difference in the Feynman rules for 4-point vertices such as VVVG. 
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and the unphysical degrees of freedom are canceled by ghost loops. The UV poles are removed by 
MS-renormalization of a s , while the IR singularities cancel against those from the real emission 
contribution (see Eq. ( fT2] ) below). 







mm" 



Figure 2: Examples of NLO QCD virtual Feynman diagrams for pp 
represent ghost particles. 




G+jet. The dotted loops 



Two independent calculations are performed for the virtual corrections. The first calculation 



is based on the Mathematica package FeynCalc [25]. Because of the Lorentz indices of the spin- 
2 graviton, we encounter high-rank tensor integrals, such as rank-5 4-point functions. Special 
care is taken to reduce those to one-loop scalar integrals by an independent Mathematica code, 
following the prescription of Ref. [^6| . 

In the second calculation, the one- loop diagrams are generated with QGRAF [27] and then 
projected onto helicity components and amplitude coefficients with FORM pq |. Details of the 
projection are given in the appendix and can also be found in Ref. [29], where the same method 



has been used. The tensor reduction is performed according to the GOLEM [30, ^] reduction 
algorithm, supplemented with additional tensor reduction routines for rank N + 1 iV-point 
tensor integrals with N < 3. To calculate numerical results we employed the OmniComp- 
Dvegas package which facilitates parallelised adaptive Monte Carlo integration and was 
developed in the context of Ref. |J3]. 



2.2 Real-emission contributions 

The real-emission contribution comprises the radiation of a real gluon or a massless (anti-)quark. 



Soft and collinear singularities are isolated using dipole subtraction [21]. Collinear emission from 
initial state partons is factorized into the parton distribution functions defined in the MS-scheme. 
The remaining IR and IR/collinear singularities, which cancel those of the virtual corrections, 
read in the notation of ||2l|] 



99 1 
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(2.2) 



(2.3) 
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where s, t, and u are the Mandelstam variables, /3o = (HCa — 4n/3jj)/3 with n/ = 5, and 
(/(e)) can be obtained by crossing from (/(e)) ^. 

The real emission matrix elements and subtraction terms are generated with MadGraph with 



Spin-2 particles [22] and MadDipole [p3f|, respectively, and are implemented in a parton- level 



Monte Carlo program. 
3. Numerical Results 

In this section we present NLO cross sections for pp/pp — > jet + G at the Tevatron (y/S = 1.96 
TeV) and the LHC (y/S = 14 TeV). Before we proceed with the numerical results, we note that 
the interaction of the KK gravitons with SM matter is described by an effective theory ||, [7j , 
which is valid only for scattering energies y/§ smaller than the fundamental scale Ms. While 
hadron collider cross sections in principle involve partonic scatterings with energies \fs up to 
the collider energy y/S, the rapid decrease of the parton luminosities at large y/§ suppresses the 
high-energy region and allows for a cross section prediction that is not very sensitive to the UV 
completion of the effective theory. We will return to this issue at the end of the section and 
provide quantitative estimates of the UV sensitivity of the theoretical prediction. 

The numerical results presented below are obtained with a s and the parton distribution 
functions defined in the MS-scheme, with five active flavours. Throughout our calculation, we 



employ the 2008 MSTW LO(NLO) PDF |3J at LO(NLO), with the corresponding value for 



the strong coupling a s . Our default choice for the renormalization and factorization scale is the 
transverse momentum of the graviton, /x = PP 



r ■ 



To suppress SM backgrounds in the LHC graviton searches [11|, we require 



P^ iss > 500 GeV . (3.1) 

The jets are defined by the hr algorithm, with the resolution parameter set to D = 0.6, and are 
required to satisfy \rjj\ < 4.5 and Pj, > 50 GeV. 

For the Tevatron predictions we use the same settings as in the recent CDF study ||, i.e. 

P^ iss > 120 GeV , P| > 150 GeV , and \r]j\ < 1 . (3.2) 

Also here, jets are defined by the algorithm with D = 0.7, and are required to satisfy \rjj\ < 3.6 
and P^ > 20 GeV. A second jet with Pt > 60 GeV is vetoed. We will, however, also discuss 
results without the jet-veto. 

We first focus on the scale dependence of the total cross section. For illustration, we set the 
model parameters to 5 = 4, and Ms = 5 TeV (LHC) and 1 TeV (Tevatron). Note that the cross 
section scales as a oc M g so that results for other values of Ms can be obtained by rescaling 
our predictions. Unfortunately, it is not possible to determine both Ms and 8 independently 
from just the shape of the P™ 1SS spectrum [fLl}] . To resolve these parameters would need very 
accurate measurements at different hadron collider center-of-mass energies; the ratio of graviton 
production cross sections at different center-of-mass energies depends on 5 through the kinematic 
limit on the graviton mass, while the dependence on Ms cancels. Operating the LHC at 7 TeV 
and 14 TeV center-of-mass energies may offer such an opportunity. We impose the kinematical 
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cuts listed in Eqs. ( |3.lD and (|3.2|) for the LHC and the Tevatron, respectively. The LO and 
NLO results are shown in Fig. | as a function of the renormalization and factorization scales 
varied around the central scale \i = P^ . We observe that the scale dependence of the NLO 
cross section is significantly smaller than that of the LO cross section, both at the LHC and at 
the Tevatron: changing n in the range between Pp /2 and 2Pjt, the LO cross section varies by 
~ 30%, while the scale uncertainty at NLO is less than ~ 10%. We have also varied both scales 
independently and find that in all cases the NLO uncertainty is less than approximately 10%. 
At the LHC, the iT-factor, K = (Tnlo/^lO) is sizeable and positive at the central scale /x = Pjt, 
increasing the LO cross section prediction by about 20%. At the Tevatron, the QCD corrections 
are mild near \i = Pp with K ~ 1, but are essential to reduce the theoretical uncertainty. 

The experimental analyses at the LHC and the Tevatron rely on the P™ 1SS and Pp* dis- 
tributions, respectively. In Figs. ||] and || we show the scale dependence of these distributions, 
for different choices of the number of extra dimensions 5 = 2,4, 6. Current Tevatron limits 
appear to exclude Ms = 1 TeV for 5 < 4. However, these analyses are based on leading- 
order predictions and should be refined using the NLO results presented here. We thus include 
numerical predictions for Mg = 1 TeV and 5 = 2, 4, 6 below. We also show the NLO QCD 
predictions for the main background pp — > Z(—> vu)+]et obtained with MCFM |J5|. [Note that 
establishing an excess in graviton plus jet production at hadron colliders requires an excel- 
lent experimental understanding of the SM background. To precisely estimate the dominant 
pp — > Z(—> i^)+jet background process one can rely on a calibration sample of the related 
process pp — > Z(—¥ e + e~ / [i + jet || [l(], 11, [DJ. Furthermore, the signal to background 



ratio can be improved by increasing the P™ lss -cut.] The bands show the uncertainty of the LO 
and NLO predictions when varying the renormalization and factorization scales in the range 
-Pjf/2 < [i < 2 Pp. The reduction of the scale uncertainty at NLO is evident. Figs. ^|and [| also 
display the Pt dependence of the iT-factors, defined as K(Pt) = (dCTNLo/dPr)/(d<7Lo/dPr)- 
The K factors are sizeable at the LHC (Fig. |J), as noted before, increasing with decreasing 8. 
Furthermore, the iT-factors depend on the kinematics and increase with increasing P™ 1SS . At 
the Tevatron, the ^-factors are in general near or below one and only mildly depend on the jet 
transverse momentum, see Fig. |5[ 

We have also investigated the uncertainty of the NLO cross section prediction due to the 



parton distribution function. Using the MSTW error PDFs [34], we find an uncertainty of less 
than approximately 15% for graviton production at the LHC, even for large P™ 1SS > 1 TeV. At 
the Tevatron, the uncertainty is even smaller and approximately 5%. 

Let us now examine the contribution of the real emission cross section with two hard jets. 
In Fig. [6| we show the ratio of the cross section where we require two hard jets with P^ > P^ ut 
and the inclusive cross section, as a function of P™ 1SS and the leading P^ et at the LHC and the 
Tevatron, respectively. Results are presented for the minimum jet- Pt set to P£ ut = 150 GeV 
and 250 GeV for the LHC, and P T ut = 60 GeV for the Tevatron. Moreover, in Fig. §, we show 
results with an alternative choice of settings (labeled B in the plot) defined as 



Pf = min(P T ) and a s = y a s (Pj}) a s (Pj?) , (3-3) 



which was used in Ref. [16] for the real emission contribution with two hard jets. We observe 
that the fraction of events containing two jets with P T > 250 GeV is 20-40% at the LHC for 
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P™ 1SS above 1 TeV, depending on the choice of input parameters for the 2-jet contribution. 
Between 40% and 70% of the events with P^ iss > 1 TeV contain two jets with P d T > 150 GeV. 
We note that even for the scale choice B, as given in Eq. (|3.3| ), the fraction of di-jet events is 
smaller than estimated in Ref. This is due to the denominator, i.e. the larger inclusive cross 
section as predicted at NLO. Nevertheless, also with our new estimate we expect a large fraction 
of high P™ 1SS events with two or more hard jets at the LHC. While the quantitative estimates 



given in Ref. [16] are thus changed due to the impact of the NLO corrections, the qualitative 
conclusions remain valid. At the Tevatron, the contribution of 2-jet events with Pj, > 60 GeV 
is moderate and does not exceed 20%. The difference between the results with the two different 
scale settings, which represent part of the uncertainty for the (tree-level) 2-jet cross section, 
increases with increasing Pt, as the difference between our default choice of scale and a s and 
the alternative choice ( |3.3| ) becomes larger. 

As mentioned at the beginning of this section, the results of the effective field theory calcu- 
lation are valid only as long as the scales involved in the hard scattering process do not exceed 
the fundamental scale Mg. To quantify the sensitivity of our prediction to the unknown UV 
completion of the theory, we compare our NLO results with those involving a truncation scheme 
which sets the cross section to zero if Qtruncation > Ms- In the numerical results presented below, 
the truncation parameter Qtruncation is taken to be the invariant mass of the missing momentum 
and observable jet(s), 

Qtruncation = | J P G + J P jet(s) |. (3.4) 

This definition is equal to Qtruncation = at LO, but takes into account that the effective 
partonic energy of the scattering process can be reduced by collinear initial state radiation 
at NLO and thus provides an IR-safe definition of the truncation parameter. NLO results 
for the transverse momentum distributions with and without the hard truncation scheme for 
5 = 2, 4, 6 at the LHC and the Tevatron are shown in Figs. [/] and [8|, respectively. As expected, 
the differences between the two calculations increase with increasing 5, as the average graviton 
mass is shifted to larger values. Also, the differences become larger as P™ 1SS (P^ *) increase. For 
example, for the LHC, at P™ ms = 1250 GeV, the uncertainties for 5 = 2,4,6 are about 5%, 
20%, and 50%, respectively, while for the Tevatron, at P£ = 250 GeV, the uncertainties for 
5 = 2,4, 6 are about 2%, 10%, and 25%. Note that the results with the hard truncation do not 
obey the simple scaling a oc M g 

Finally, we comment on the prospects for graviton searches during the initial phase of the 
LHC operating at 7 TeV. Even at half the nominal center-of-mass energy, the LHC will be 



able to extend the sensitivity of current Tevatron searches to larger values of Ms 13, 14]. For 
illustration we show in Fig. ^ the NLO P™ 1SS distribution for Mg = 2 TeV and (5 = 2,4, 6, 
together with the dominant SM background. We find sizeable signal rates, in particular for 
5 = 2, which exceed the background for P™ 1SS > 250 GeV. We re-emphasize that comparing 
mono-jet signatures at 7 TeV with results obtained at higher center-of-mass energies during a 
later stage of LHC operation may allow to disentangle the fundamental parameters Ms and 5. 
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4. Summary 

We have presented the first calculation of the NLO QCD corrections to Kaluza-Klein graviton 
plus jet hadro-production in models with large extra dimensions. The calculation has been set 
up as a fully-flexible parton- level Monte Carlo program 2 , and results have been presented for 
cross sections and distributions at the Tevatron and at the LHC. 

The QCD corrections stabilize the theoretical prediction and significantly reduce the scale 
uncertainty to a level of approximately 10%. Near the central scale, [i = P!f , the QCD correc- 
tions increase the cross section at the LHC by 30-50%, depending on the kinematical region and 
the choice of model parameters. At the Tevatron, the QCD corrections are modest and negative 
near /x = PS and do not strongly depend on the kinematics. A significant contribution of di-jet 
events is expected at the LHC, where 20-40% of signal events at P™ 1SS > 1 TeV contain two jets 
with Pj, > 250 GeV. At the Tevatron, on the other hand, the contribution of 2-jet events with 
Pj, > 60 GeV is moderate and does not exceed 20%. The theoretical uncertainty of the cross 
section prediction due to the parton distribution functions is mild, with approximately 15% and 
5% at the LHC and the Tevatron, respectively. 

We have also studied the uncertainties arising from the UV completion of the theory by com- 
paring our default NLO results with those involving a hard truncation scheme. The differences 
between the two calculations are small for 5 = 2 but can reach up to 50% for 5 = 6 and large 
Pt- Reducing these uncertainties requires to go beyond the effective field theory approximation 
of Eq. ( |1.2p , which is beyond the scope of the present paper. Ignoring form factor effects for 
the graviton couplings to gluons and quarks and ignoring Kaluza-Klein excitations in the loops, 
defines one particular phenomenological model. For this model, our calculation quantifies the 
size of QCD corrections, and these results may then be taken as an indication of what to expect 
of QCD corrections in more complete models of the UV physics. 



Acknowledgments 

This work is supported in part by the Deutsche Forschungsgemeinschaft under SFB/TR-9 "Com- 
putergestiitzte Theoretische Teilchenphysik" , the Helmholtz Alliance "Physics at the Terascale" , 
and the European Community's Marie-Curie Research Training Network under contract MRTN- 
CT-2006-035505 "Tools and Precision Calculations for Physics Discoveries at Colliders". SK is 
grateful to Thomas Binoth and Nikolas Kauer for valuable discussions. SK and MK would like 
to thank Arnd Meyer for discussions on experimental aspects of mono-jet searches. 



A. Amplitude calculation 

The amplitude for the partonic processes gg — > Gg and qq — > Gg can be expressed as 

M(g( Pl , X 1 )g(p 2 , A 2 ) -»■ g(p 3 , X 3 )G(p 4 , A 4 )) = M^^e^e^ f abc 
■M(9(pi,Ai)g(p2,A 2 ) ^<7(p 3 ,A 3 )G(p4,A4)) = ip^\M^ v \p^) e^e**!*, 



The fortran code is available upon request from karg@physik.rwth-aachen.de. 
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where \p i ) is the Weyl spinor for a massless particle with momentum pi. Since we consider 
all quarks to be massless, the helicity of the quark line is conserved, and therefore we have 
Ai = A2 G [+, — ]. (Note that the physical helicity of the anti-quark is given by — A2.) 

Applying spinor helicity methods |36|], the polarization vector for a massless spin-1 boson 
for helicity ± is given by 

4(P.r) = ± { -^4, (A.2) 



where the vector r in ( [A. 2 ) denotes an arbitrary reference vector (with r 2 = 0). 

The polarization tensor e^, of the graviton, a massive spin-2 vector boson, can be con- 
structed from the polarization vectors of massive spin-1 bosons, en' , as follows: 3 



e ++ = e + e + 
e+ = — (e + e° + e°e+ 

y/2 i-' . 



1 



eL = ^{e+e- + e-et-2ey v ) (A.3) 



e~ = — (e-e° + eV) 

In order to extend the spinor-formalism to massive gauge bosons, it is useful to decompose 
the graviton momentum into two light-like vectors: 

771 2 

P4 = qA + a r, with p\ = m 2 , q\ = = r 2 , a = , (A. 4) 

2 p^ ■ r 

where the arbitrary reference momentum r can be taken from the list of available light-like ex- 
ternal momenta. The expressions for the three polarization vectors e ±,Q can now be constructed 
from the two light-like vectors p^ and r and read 

£ ±( m m) - ± ( rT \^\lh 

V2(rT\p±) (A5) 

1 

eu(P4,m) = — (<?£ -ar 11 ) . 
p m 

Of course, individual helicity amplitudes are no longer independent of the choice of the 
reference momentum of the graviton. However, we are only interested in the spin sum, which is 
independent of the reference momentum. 

By a suitable choice of the reference vectors, we can assemble the individual spinor products 
and write them as a trace times a global spinorial factor. The projector for the A1A2A3 = + + + 
helicity combination, for example, reads: 



.+ r + .+ _ (3 |mi 1 1 ) (1 l/fgjj ) (2 |/x 3 |3 ) = tr [(1 -7 5 )3/xil/x 2 2/i 3 j , A Q , 

e w e M2 e M3 ^2(31) V2 (12) V2(23) 4^2 (31) (12) (23) ' , "" > ' 



Note that Eq. (A. 2) follows the convention e*'* = +e„ . An additional sign is sometimes included in the 
literature, which would alter the sign of the third term in e°„ in Eq. (A.3) 
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and the projector for the A4 = 2+ helicity is given by 



++= (l-H4-)(4-|^ 2 |l-)(l-H4-)(4-|^ 2 |l-) _ tr [(1-75)1^42^42] 
V " V2 (14) (42) [21] ^2 (14) (42) [21] 4 (14) 2 (42) 2 [21] 2 ' 1 ' ' 

with the spinor inner products (ij) = (p^\p1f), [ij] = (pf\pj)- 

Note that the 40 (20) helicities for gg — > gG (qq — > gG) are related to each other by discrete 
symmetries, like parity, Bose symmetry or invariance under charge conjugation. The symmetries 
therefore allow for a cross check of the results or can be used to reduce the algebraical work by 
calculating only a generic set of helicity amplitudes. In this way, we could perform the spin sum 
by computing only 1 (2) helicity amplitude(s) for gg — > gG (qq — > gG). 

The helicity amplitudes for gg — > gG and qq —> gG summed over the polarization of the 
graviton are given by 

£ | - M ^' 2 = £ and E '^cl 2 = ^ (4*1 + sm%) (A.8) 
A4 A4 

where u = u — uiq etc. Applying Bose and parity transformations on the above expressions 
gives us all helicity amplitudes. 

Summing over final colours and averaging over initial helicities and colours, the squared LO 
matrix elements are given by 



\M LO \\gg -> gG) 1 ' 



32 Mp V" stu 



p 



\M LO \ 2 (qq^ gG) = -^r[(4tu + sm 2 G ) 

9 Mp V stu 



(A- 



in agreement with ||, but in a form where the symmetries are more obvious. The process 
qg — > qG is related by crossing to qq — > gG. 
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Figure 3: Scale variation for the integrated cross section at LHC and Tevatron, for a common scale 
fi = /j, r = fif. Selection cuts are described in the text. 
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Figure 4: P™ lss distribution for the graviton signal at the LHC with scale uncertainty bands (0.5 Pj! < 
fi < 2Py? ). Also given is the NLO distribution for the dominant Z — > vv background. The lower part of 
the plot shows K(P T ) = (dcr N Lo/dPT)/(dCT L o/dPT) for 5 = 2,4,6 (top down). 
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Figure 5: Same as Fig. [| but for the Pt distribution of the leading jet at the Tevatron. 
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Figure 7: Effect of truncation of the partonic cross section above (^truncation = Ms at the LHC. See text 
for details. 



10 



0.1 



dcr r fb 



dP ^t LGeVJ 



Tevatron 
M s = 1 TeV 



NLO 
NLO trunc 




150 200 



250 300 
P^ et [GeV] 



350 400 
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